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We propose a systematic scheme to reach the properties of two-dimensional (2D) statistical and quantum 
systems by studying the effective (l-l-l)-dimensional theory that is constructed from the tensor network repre¬ 
sentation. On on hand, we discover that the degeneracy of the 2D system can be determined by the purity of the 
boundary thermal state, which is the density operator of the effective theory at zero (effective) temperature. On 
the other hand, we find that the gapped (or critical) 2D system leads to a gapped (or critical) effective (l-l-l)- 
dimensional theory whose criticality can be accessed by the entanglement entropy S of its ground state dubbed 
as boundary pure state. We also uncover that for the critical systems, S obeys the same logarithmic law as that 
found in the critical ID quantum chains, which reads S = (kcI6) logj D + const., with c the central charge and 
K a constant related to the scaling property of the correlation length f as ^ ~ D". Such a scaling law presents an 
efficient way to characterize the critical universality class of the original 2D systems. An important implication 
of our work is that many well-established theories for ID quantum chains become available for studying 2D 
systems with the help of the proposed lower dimensional correspondence. 

PACS numbers: 71.27.+a, 74.40.Kb, 03.65.Ud 


Introduction .— Lots of efforts have been made on explor¬ 
ing novel properties of quantum many-body systems in two 
dimensions (2D), where the fruitful geometries of the 2D lat¬ 
tices and the strong competition between the quantum fluctua¬ 
tion and magnetic ordering provide a fertile ground for various 
exotic phenomena. For instance, the 2D frustrated Heisenberg 
models such as the kagome antiferromagnet I!]] were shown to 
be good candidates for realizing quantum spin liquids (QSL’s) 
0. It is of basic scientific interest to study the elusive prop¬ 
erties of the QSL’s such as topological orders 0, fractional 
excitations |3l and criticality 0. 

However, many important issues are hard to access in 2D. 
Because the bipartition of a 2D system suffers much more 
complexity than that in one dimension (ID), quantities which 
involves the subregion size scaling, e.g. the central charge 
which characterizes the critical universality class fill and the 
topological entanglement entropy which can be used to iden¬ 
tify topological orders fi], are extremely difficult to reach. A 
usual way to make a compromise is to calculate on a finite 
torus lUilsl], which suffers the finite size effect. Meanwhile, 
classical simulations are inefficient on critical systems even in 
IDllI^. which makes the criticality in 2D more elusive. 


The boundary theories that are recentl y pr oposed based on 
the tensor network (TN) representation ll iMoll present us a 
new clue to solve problems that cannot be easily handled with 
conventional methods. The tensor network state (TNS, also 
termed as projected entangled pair state) iSilll gives a nat¬ 
ural representation for a quantum state, where the local tensors 
act as the projectors that map the physical degrees of freedom 
to virtual ones carrying the entanglement. The inner product 
of a TNS and its copy {<fTNs\^TNs) maps a quantum many- 
body state to its classical correspondence in the same dimen¬ 
sion. The resulting classical model can be further mapped 
onto an effective ID quantum model. Some theories based 


on this scheme have been proposed to tackle difficult prob¬ 
lems, e.g. identifying the topological orders 11211 and simulat¬ 
ing low-energy excitations 11311 of 2D quantum systems, etc. 

In this work, we propose a correspondence from the TN 
representation of a statistical system or a quantum state in 2D 
to an effective (l-i-l)-dimensional theory ff. Different from 
the existing theories that are based on the reduced density ma¬ 
trix, ff is constructed straightforwardly by the transfer matrix 
of the TN. We suggest that the properties of the original 2D 
system can be efficiently reached by investigating this lower 
dimensional correspondence ff. 

First, the boundary thermal state (BTS) is introduced as the 
thermal state of ff at zero effective temperature, which is ob¬ 
tained by the linearized tensor renormalization group (LTRG) 
method 0. Note the LTRG was proposed to obtain the ther¬ 
modynamic properties of ID quantum systems. We find the 
purity of the BTS a robust quantity to determine the degener¬ 
acy of the 2D system, i.e. the BTS is pure if the corresponding 
parent Hamiltonian llTIl is non-degenerate, and is mixed if the 
parent Hamiltonian is degenerate. 

Second, the boundary pure state (BPS) defined as the 
ground state of the (1 -H1 )-dimension theory can be obtained by 
the infinite time-evolving block decimation (iTEBD) 01811191] 
method. Note that the iTEBD was suggested for calculating 
the ground states of ID quantum chains. We discover that the 
entanglement entropy {S) of the BPS can be utilized to de¬ 
tect the criticality of the system. Eor gapped systems, the BPS 
is found to bear only a finite S, while for critical systems, 
S increases with the dimension D of the BPS, obeying the 
same logarithmic law as that found in the critical ID quantum 
chains llOll. which reads 


S — — log 2 D + const., 
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with c the central charge (0 and k a constant related to the 
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scaling of the correlation length ^ defined as ^ ~ D* iI3]. 
This gives an efficient way to determine the critical universal¬ 
ity class of the 2D statistical and quantum systems through 
its (l-i-l)-dimensional correspondence, while the difficulties 
in scaling against the subregion length in 2D are avoided. 

To examine our proposal, we perform exact deductions 
and numerical simulations on several celebrated examples in¬ 
cluding the 2D Ising model at and away from the critical 
temperature, the Greenberger-Horne-Zeilinger (GHZ) state 
(gapped and non-topological) 111 ], the Z 2 state (gapped and 
topological) 13], the nearest-neighbor resonating valence bond 
(NNRVB) state on kagome (gapped and topological) and hon¬ 
eycomb (critical and topological) lattices 12211 . 

Correspondence from a 2D TN to a (1+1 )-dimensional the¬ 
ory .— A planar TN is defined by the contraction of tensors 


Z = 


= 'Cf,). 


( 2 ) 


where . is the local tensors, j runs over all tensors and 
tTr is the trace of all common bonds. Eq. yji plays a central 
role in the TN scheme since the calculations of many physi¬ 
cal quantities (e.g. the partition function of the 2D classical 
models 123 or the inner product of two quantum states jfSj) 
is equivalent to computing such a TN contraction. Here we 
take the square TN as an example [Fig. [T|(a)]. 

Let us first explain how to apply the LTRG and iTEBD 
methods to calculating Eq. (|2]i. Without losing the general¬ 
ity, we presume the TN satisfies the translational invariance, 
i.e. . - Tuidr- We start the contraction with the original 

lljljdjrj 

local tensor E® ^ = TuUr- For the fth step of the contraction, 
the effective tensor is updated by contracting one original 
tensor TuUr to as = Y^du’ TiddrP^i7^l,,,ddu' with 

I" - (I, I') and r" = (r, r'), as shown in Fig. [T](a). Such a con¬ 
traction scheme is equivalent to applying the LTRG method to 
contract a 2D planar TN instead of performing the imaginary 
time evolution of the ID quantum systems in its original pro¬ 
posal 116 f The TN is contracted to an matrix product operator 
(MPO) 12411 formed by and the entanglement spectrum T® 
residing on the virtual bonds. The dimensions of I" and r" of 
P are bounded in the same way as that the LTRG bounds the 
dimensions of the MPO, i.e. the truncations are implemented 
according to the entanglement spectrum after canonical- 
izing the MPO. The local canonical conditions of an MPO, 
which is similar to the canonicalization of an matrix product 
state (MPS) Il3 . are shown in Fig. [I](b). Specifically speak¬ 
ing, an MPO is canonical when the spectrum on any bond is 
the singular value spectrum between the left/right subregions. 

At the fth step, represents efficiently a tensor stripe 
containing t + 1 original tensors [dash squares in Fig. [1] 
(a)], and the MPO represents the TN with f H- 1 layers. As 
f —> 00 , the MPO converges to the fixed point, i.e. PuVd'r" - 
Ydii' TuidrPu’i’d’r’Sdu’- The fixed point MPO dubbed as the BTS 
represents effectively the whole TN. 

Such a contraction of a planar TN can also be done accu¬ 
rately using the MPS with the iTEBD algorithm jJA, 13 ■ Af¬ 



FIG. 1: (Color online) (a) The LTRG scheme on a square TN formed 
hy the local tensor T and its copies. The contraction of the TN is 
performed linearly with the help of an MPO formed hy local tensors 
P and spectrum A^. At the fth renormalization step, each local tensor 
effectively represents a tensor strip (the dashed circles) containing 
(f -I- 1) original tensors. The transfer matrix of the TN is defined 
as the tensor stripe along the horizontal direction, (h) The left and 
right canonical conditions for the MPO. 


ter making a sufficiently large number of contractions, one 
obtains the converged MPS with its local tensor and entangle¬ 
ment spectrum denoted by A and . We dub such an MPS 
that is pure as the BPS of the TN. 

The contractions of the TN through the LTRG and iTEBD 
are actually equivalent to the imaginary time evolution of a 
(l-i-l)-dimensional theory. To see this, we introduce the trans¬ 
fer operator of the TN defined by the infinite tensor stripe 
along the direction perpendicular to the contraction direction 
[the yellow shadow in Fig. [I](a)]. The Hamiltonian p{ of 
the (iH-l)-dimensional theory then can defined as e = p™. 
Here, we require the TN to have spatial inversion symmetry 
so that Pi is guaranteed to be hermitian. 

In the LTRG, we readily have the MPO at the fth step as 
P^MPO ~ e with “TC the effective inverse temperature sat¬ 
isfying 'K - t + \. In the iTEBD, contracting the fth layer 
of the tensors in the TN to the MPS is actually equivalent to 
evolving the MPS along the effective imaginary time direc¬ 
tion as \MPS), = e-'^\MPS),-i = e-'^'^\MPS)o with \MPS)o 
the initial MPS. The BPS defined as the fixed point MPS is 
the ground state of Pd, fulfilling e^'^\BPS) - C\BPS) (C is a 
constant) which is just the fixed point condition. 

Purity of the boundary thermal state .— Now we show that 
the degeneracy of the ground state(s) of Pd can be detected by 
the purity of the BTS, and suggest that such a degeneracy can 
identify the degeneracy of the original 2D system. When 
is non-degenerate, the BPS \ fBPs) is simply the ground of Pd, 
and the BTS pbts gives a pure state that is the outer product of 
\il/BPs) audits copy, say p^j-j = \ fBPs){fBPs\- Consequently, 
the entanglement spectra and of the BTS and BPS have 

a simple outer-product relation A'^ - A‘^ ( 8 > A‘^. 

When Pd is degenerate, the BTS is no longer pure, i.e. 
it cannot be decomposed into an outer product form. The 
BTS can be formally written in a mixed thermal state pBPs = 
with 10 ') the ith degenerate eigenvector, 77 the 
thermal probability distribution and x the degeneracy. Be¬ 
cause the degenerate eigenvectors of p^*^ correspond to a 
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same eigenvalue, the BTS given by LTRG (which is essen¬ 
tially a power method) is expected to be the maximal mix¬ 
ture of the degenerate states, i.e. rji - 1/x- Meanwhile, it 
is widely accepted that the BPS \iJ/bps) by iTEBD favors the 
minimally entangled state among all combinations of the de¬ 
generate eigenvectors. In this case, the entanglement spectra 
of the BTS and BPS do not have a simple outer product rela¬ 
tion. Thus we have a way to identify the purity of the BTS by 
monitoring the spectrum difference s - \A^ — |. 

We observe that the degeneracy of "N is in accordance with 
the degeneracy of the original 2D system. We shall remark 
that by the degeneracy of a quantum state, we mean the de¬ 


generacy of its parent Hamiltonian Ill7l| (n ote a statistical sys¬ 


tem can be mapped to a quantum state Il25[] l. Interestingly, the 
BTS gives a mixed state whether the parent Hamiltonian is 
non-topological (e.g. the GHZ state) or topological (e.g. the 
Z 2 topological state). Thus the purity of the BTS is robust to 
detect the degeneracy of the system under interest. 

Entanglement scaling of the boundary pure state. — Many 
remarkable works have been done with the entanglement of 
the MPS to study the properties of the ground states of ID 
quantum models; the subregion length-dependence of the en¬ 
tanglement of a critical state of a spin chain is analogous to 
that of the entropy in conformal field theories ifioll : the en¬ 
tanglement exhibits the correlation properties and the critical 
exponents lE^ of the system as well as the topological prop- 
erties, such as the symmetry-protected topological orders 11271] 
and edge excitations El. Through our correspondence, these 
achievements become available for investigating 2D systems. 

We suggest that the original 2D system shares the same crit¬ 
icality with its lower dimensional correspondence 'H. Speak¬ 
ing in detail, for a gapped 2D system, the TY is gapped where 
the BPS (the ground state of TY) possesses only finite entan¬ 
glement entropy from the known conclusions in ID quantum 
chains. For a critical 2D system, the TY is critical, where the 
entanglement entropy defined as S (d'^) = - XfAff log 2 (Aff 
of the BPS obeys a logarithmic scaling law given by Eq. O- 

To study the criticality of (infinite) strongly correlated sys¬ 
tems directly in 2D is a extremely hard task both theoretically 
and numerically: (a) the conventional subregion scaling for 
the central charge S (A^) = (c/3) log 2 L H- H is difficult to 
implement in 2D since the bipartition is essentially different 
from that in ID; (b) any finite-entanglement approximation of 
a critical system gives an effective gapped model whose cor¬ 
relation length is finite. We speculate that the central charge 
of TY which can be efficiently obtained ifl^ can characterize 
the critical universality class of the original 2D system. Thus 
the difficulties encountered in 2D can be avoided. 

Two exactly contractible states. — We give the exact deduc¬ 
tion on two non-trivial states, the GHZ 111] and Z 2 topolog¬ 
ical state 13] to examine our proposal on detecting degener¬ 
acy. The GHZ state is a topological-trivial state whose parent 
Hamiltonian has two degenerate ground states. The Z 2 state is 
the exact ground state of the toric-code model which is gapped 
and possesses nontrivial topological degeneracy ll29l] . 

The local tensor of the TN in Eq. (|2|i of the GHZ 


and Z 2 states can be written in a unified form TuUr - 

with A the rank 

spectrum. We have Aq = Ai = 1/ V2 and t/ a (2 x 2) 
unitary matrix S. I is the super-diagonal tensor satisfying 
Ifipi- - 1 if /i = //I - ■ ■ ■ or - 0 otherwise. The matrix 
U is the rotation matrix satisfying Uo.o = -t/ 1,1 = cos 6 and 
Uo.i = t/o.i = sin 0 with 6 the parameter. By taking 0 = 0 and 
7r/4, U becomes the identical and Fourier matrix, which gives 
the TN for the GHZ and Z 2 states, respectively. 

The BTS of the GHZ and Z 2 TN can also be ex¬ 
pressed in a unified form with the local tensor = 

2^11^12 IpAnrUufii Udf 42 and the spectrum A^ - A Q- It is easy 
to check that the BTS is canonical. Obviously, such a BTS is 
mixed and possesses a finite bond dimension D - 2 with the 
spectrum A^ - A^ - Ij \/2. For the BPS, the MPS converges 
to an unstable fixed point Aut - 2/2 Iiih'Uufi and A^ = A only 
by choosing = / as the initial MPS. We have the entan¬ 
glement entropy S (A^) = log 2 2=1 for both the GHZ and Z 2 
states. Otherwise the fixed point flows to the stable fixed point 
with zero entanglement of the MPS. 

It is not surprising to have only D = 2 of the BTS and 
BPS for the GHZ state, as it is just the superposition of two 
classical ferromagnetic states. In contrast, for the Z 2 state, 
the quantum entanglement entropy (of the physical degrees 
of freedom) obeys the area law and increases unboundedly as 
S - aL - j with L the boundary length of the subsystem and 
y the topological entanglement entropy 13 . Consequently, a 
large bond dimension ~ 2'^ ~ 2“^ is expected to capture the 
quantum entanglement directly in 2D. Amazingly, the ground 
state of its lower dimensional correspondence only has a small 
dimension D - 2 with a finite entanglement entropy. 

We show that the absence of the physical degeneracy leads 
to a pure BTS. If we destroy the degeneracy of the ground 
states of the parent Hamiltonian, we have, equivalently, a shift 
of A as A 2 = eAi (0 < e < 1). Then for the fth step of 
renormalization, the canonical spectrum of the MPO satisfies 
A^ = As f ^ 00 , dj vanishes and the BTS becomes 

a pure state. This picture also holds for the Z 2 state. 

The 2D Ising model. — The partition function of the 2D 
antiferromagnetic (AF) Ising model on square lattice can be 
written as the contraction of a TN as Eq. (|2]). The local tensor 
reads Tuw - exp{-[s„s; - 1 - sisj + sjSr + SySu + h^isu - si + sj- 
■S)')]/T}, where s;, Sd and s,- are the four spins in a plaquette, 
hs is the staggered magnetic field and T is the temperature. 

We study the purity of the BTS at different temperatures by 
numerically calculating e = \A'^ - A^ ® | shown in Fig. 

(a). For T < Tc with T^ the critical temperature, we calculated 
e with (or without) a small staggered field hs = lO^'* where 
the Z 2 symmetry breaking is triggered (or conserved) in the 
LTRG scheme. We have d® = 1, d®^j = 0 and d'^ = A^ with 
hs - 10 which means the BTS is pure (s = 0). For hs - 
0, the system contains two degenerate states (Z 2 symmetry 
is conserved) and the BTS flows to a GHZ-like mixed state 
for T < T" [dO = df = 1/V 2 , A% < 0(10-^), ^ 1, 

4t < (9(10 ^)] and e is non-zero. 
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FIG. 2: (Color online) (a) The temperature-dependence of e = \A^ - 
(g> I of the 2D Ising model with and without a small staggered 
magnetic filed /i,. The inset of (a) gives the entanglement entropy S 
of the BTS versus its bond dimension D at T = 2, 2.5. (b) At T = T^., 
the correlation length ^ scales with D as ^ ~ ly with k = 2.021. The 
inset of (b) presents that 5 satisfies S = ki logj D + const., with ki = 
0.166. According to Eq. Q, we have the central charge c = 0.493. 


FIG. 3: (Color online) The scaling of (a) the entanglement entropy 
5, (b) the correlation length ^ of the honeycomb NNRVB state and 
(c) the entanglement entropy 5 of the kagome NNRVB state against 
the bond dimension D of the BPS. For the honeycomb NNRVB state 
which is critical, S obey a logarithmic scaling law S = k\ log, D + 
const, with ki = 0.241 and ^ satisfies ^ ~ ly with k = 1.475, which 
gives the central charge c = 0.98 ~ 1 from Eq. ([TJ. Eor the kagome 
NNRVB state, S saturates with D showing such a state is gapped. 


For T > T‘^, the MPO flows to a trivial disordered state 
with or without a staggered field, and the BTS gives a pure 
state with vanishing e. Consequently, the separating point of 
the two curves of e gives the critical temperature accurately, 
where we have Tc - 2.27, in comparison to the exact criti¬ 
cal temperature = 2/ln(l -i- V2) ^ 2.269. One can also 
see that at both sides of the critical temperature, the entropy 
S (/l'^) saturates to a finite value when the bond dimension D 
increases, as shown in the inset of Fig. |2](a). 

At T = Tc, s is non-zero indicating the BTS is a mixed 
state. We remark that the purity of the BTS at the critical point 
is not robust as an MPO with finite D cannot accurately give 
true BTS of the TN. Fig. |2](b) gives the correlation length ^ 
versus D, which satisfies ^ ~ D'^ with k - 2.021 indicating the 
(l-i-l)-dimensional theory 'K is critical in. Such a criticality 
is also supported by the logarithmic scaling law of S of the 
BPS shown in the inset of Fig. |2](b). Then we accurately have 
the central charge c - 0.493 from Eq. O, while the exact 
result is c = 1/2 from the conformal field theory fl. 

The topological RVB states .— The TNS representations of 
the NNRVB states are given in Ref. 02511 . It is known that the 
NNRVB, which is topological, is critical on bipartite lattices 
03 in but gapped on non-bipartite lattices 0220 . 

Fig. [3](a) shows that the entanglement entropy S of the BPS 
of the honeycomb RVB state increases with D in a logarithmic 
way as Eq. O with the coefficient (a:c)/6 = 0.241, and we 
have K = 1.475 with the scaling of the correlation length ^ 
shown in Eig. [3](b). Then we have the central charge of the 
NNRVB on honeycomb lattice c - 0.98 « 1. 

Eor the kagome RVB state which is gapped, we find s ^ 
0.1, indicating that the BTS is mixed. In Eig. |3](c), we observe 
S saturates at about S ^0.19, evidencing that our theory can 
identify whether a 2D quantum state is gapped or gapless by 
the criticality of the (iH-l)-dimensional correspondence. We 
do not see the degeneracy in the entanglement spectrum of the 
BPS, which we have however seen in the Z 2 state. The reason 


may be that the corresponding symmetry is not well protected 
during the contraction procedure 13211 . 

Conclusion .— In summary, we propose that the degener¬ 
acy and criticality of a 2D system can be efficiently identi¬ 
fied by the purity of the BTS and the entanglement entropy of 
the BPS, both of which are associated to a lower dimensional 
correspondence 'K. We suggest that the central charge of 'K 
which is reached by the finite dimension scaling can be used 
to characterize the criticality of the original 2D systems. Ex¬ 
act deductions and numerical calculations on different kinds 
of examples give strong supports to our proposal, where we 
uncover that the NNRVB state on honeycomb lattice bears the 
central charge c a; 1 . 

Our proposal could be readily extended to study other 
critical phenomena and critical quantum states such as the 
entropy-driven phase transitions i33n . critical quantum lattice- 
gas states 13411 and chiral spin liquids 1351] . as well as the vari¬ 
ational ground states obtained by the TN-based algorithms. 

This work is supported in part by the MOST of China 
(Grant No. 2012CB932900 and No. 2013CB933401), the 
Strategic Priority Research Program of the Chinese Academy 
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Supplemental Material 

I. Tensor networks of the GHZ and Z2 states 

The GHZ [[ll] state is a highly entangled state 
whieh is defined as |iAg//z) = -;) 5 (nZi I0)i + 

n^i |l)y) and has been introduced as a source state 
for quantum computations 
form of a square TNS is 


The GHZ state in 


\^GHz) - zn 4„,ia„,2a„,3£!„,4 Y\ 


( 1 ) 


|a,i| n 


where \sj) with sj = 0, 1 denotes the up or down 
eigenstate of the jth spin and {a} are the ancillary 
indices. I is the super-diagonal tensors defined as 




aia2-a„ 


1 , ai = a2 = ■ 

0, otherwise. 


a„. 


( 2 ) 


The Z 2 topological state is the ground state of 
Z 2 Hamiltonian Hz^ = -U ZiUieiegs of - 
t Zp W judges of p ^ *e Pauli op¬ 

erator yD. It is a topologically ordered state with 
long range entanglement. Its TNS representation is 

I*Az 2) = ]~| 2a„,ia„,2a„,3an,4 ]”[ ^SjOjcA^j)^ (3) 

|a,i| n i 


where Q is 


Q 


^n,l^n,2^n,3^n,4 


1. ZLi ^n,a = even. 
0, otherwise. 


(4) 


One can see that the TNS of such two states can 
be written in a unified form 


|a,i| n 


pin) 


n 


rU) 




(5) 


where is a tensor defined on the nth vertex of 
the network and is a super-diagonal tensor [Eq. 
([2])] on the jth edge. See FIG. [STIfaf. The tensor 
is actually a projector which projects the ancil¬ 
lary indices represented by a, and a', to the physical 
index sj. 

With the given TNS, a calculation of concerned 
quantity, such as Z = ((A|iA) or (O) = {il/\d\i//)jZ 
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FIG. S2: (Color online) The sketch that shows the fixed point MPO. 
For the first arrow, we apply the orthogonality UU^ = /, and for the 
second arrow, we apply the equations - ~ 


FIG. SI: (Color online) (a) The unified TNS representation of the 
GHZ and Z 2 state [Eq. ©]. (h) The sketch of the inner product TN 
Z = (i/'|(/'). (c) The TN Z can be transformed with Eq. into a 
TN formed by the local tensor F shown in (d). (e) F can be decom¬ 
posed by the orthogonal tensor decomposition with Eq. (f) The 
eigenvalue decomposition of the tensor I [Eq. 0]. 


with 0 a quantum operator, becomes the contrac¬ 
tion of a corresponding TN. For the TN of Z, the 
local tensor satisfies [FIG.[ST](c)] 

8n,l§n,2§n,3Sn,4 / J ^n,l^n,2^ii,3^n,4 ^n,l^A7,2^fz,3^/i,4 

a,b (6) 

j(n,2) j(n,3) j(nA) 

^n, 2 ^n, 28 n ,2 ^n,3^n,38n,3 ^n,4^n,4^«,4 ^ 

where (/ = 1,2,3,4) is obtained by the 

^nj^nj8nj 

eigenvalue decomposition [FIG.[Sl](f)] 


E 


j(nJ) 


ri’bj) 




E jinJ) j(n,j) 


(7) 


When taking the local tensor as the one in 
Eq. ([U) or Eq. (|3]), one readily has in 

Eq. by substitutions (or according to the tensor 
fusion algebra [0]). Meanwhile, we found that 
can be decomposed by orthogonal tensor decompo¬ 
sition dH] as 


rp(n) 

0,1,10,1,2 Un,3 0,1,4 


1 


r=0 


u, 


a„Ar, 


( 8 ) 


The matrix U for GHZ and Z 2 TN is actually the 
rotation matrix 

^ [ cos(0) since) 1 
[ sin(e) - cos(e) J 

with 9 the parameter. U is orthogonal which satis¬ 
fies 

Irr' = Yj UarUar'. ( 10 ) 


By taking = A = [1/ V2,1 / V2] in Eq. ([8]) and 
e = 0 in Eq. d^, one obtains the GHZ state with Eq. 
d5]) and the corresponding local tensor of Z with Eq. 
dH]). By taking A = [1/V2,1/V2] and 6 = 7r/4, 
one has the Z 2 topologically ordered state and the 
local tensor of Z. Comparing these two states, both 
of them are highly entangled, while the GHZ state 
was introduced for quantum teleportation and the 
Z 2 state was revealed to bear a non-trivial topolog¬ 
ical entanglement. But against expectations, such 
two states belong to the same class characterized 
by only one parameter 6. 

Note that any translational invariant TN can be 
exactly contracted when the local tenor bears the 
same form as the that of the GHZ/Z 2 TN shown 
above, which actually gives the orthogonal tensor 
decomposition dH] ; the contraction properties of the 
Z 2 state has been studied from a different perspec¬ 
tive using the Hopf algebra [0]. 

II. The fixed-point MPO of the GHZ/Z2 tensor 
networks 


EIG. |S2] illustrates the proof of the fixed point of 
the GHZ and Z 2 TN with ETRG. The local tensor 
of the MPO satisfies 


01020304 


1 

r=0 


oirIo2rUa3 r^o4 r ; 


( 11 ) 


and the entanglement spectrum is = A. Dur¬ 
ing the contraction, the adjacent matrices f/’s 
vanish according to Eq. (fTOl) . and by apply¬ 
ing E/, — Ijiyfir-- nnd E/ 2 / 2 ' 7/2/2'...//,/,'... ~ 

E /2 7/2...//,..., the resulting MPO is exactly the same 
as the one before the contraction. This proof holds 
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for both GHZ and Z 2 TNs by taking 6 = 0 and n/A 
in U [Eq. Q], 
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